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FLUTTER AND OSCILLATING AIR-FORCE CALCULATIONS FOR AN AIRFOIL IN A 

TWO-DIMENSIONAL SUPERSONIC FLOW 

By I. E. Gaerick and S. I. Rubinow 


SUMMARY 

A connected account is given of the Possio theory of non- 
stationary flow for small disturbances in a two-dimensional 
supersonic flow and of its application to the determination of the 
aerodynamic forces on an oscillating airfoil. Further applica- 
tion is made to the problem of wing flutter in the degrees of free- 
dom — torsion, bending, and aileron rotation. Numerical tables 
■for flutter calculations are provided for various values of the 
Mach number greater than unity. Results for bending-torsion 
wing flutter are shown in figures and are discussed. The static 
instabilities of divergence and aileron reversal are examined 
as is a one-degree-of-freedom case of torsional oscillatory 
instability. 

INTRODUCTION 

The problem of flutter or aerodynamic instability for 
high-speed aircraft is of considerable importance and hence 
interest is directed to the aerodynamic problem of the oscil- 
lating airfoil moving forward at high speed. Although for 
conventional aircraft the subsonic and the near-sonic or 
transonic speed ranges are still of main interest, the super- 
sonic speed range is becoming increasingly significant. 

A theoretical treatment of the oscillating airfoil of infinite 
aspect ratio moving at supersonic speed has been given by 
Possio (reference 1). This treatment is based on the theory 
of small perturbations to the main stream, thus is essentially 
an acoustic theory, and leads to linearization of the equation 
satisfied by the velocity potential. The airfoil is therefore 
assumed to be very thin, at small angle of attack, and the 
flow is assumed nonviscous, unseparated, and free from 
strong shocks. 

The small-disturbance linearized theory, being much less 
complicated than a more rigorous nonlinear theory, is to be 
regarded as an expedient which allows an initial theoretical 
solution. The theory permits the occurrence of weak (in- 
finitesimally small) shocks and thus the basic trends and 
effects of the parameters of the simplified problem can be 
indicated. The theory reduces to that of Ackeret in the 
stationary (static) case and, like it, is not expected to be 
valid too near M= 1 . In view of the restrictions and as- 
sumptions in the analysis, important modifications may be 
required in certain cases for thick finite airfoils; but even 
here the simple theory for thin wing sections may serve as a 
basis. 


In addition to Possio’s brief work, an equivalent extended 
treatment has been given by Borbely (reference 2) which 
utilizes contour integrations to carry out the solution of the 
partial differential equation for the velocity potential accord- 
ing to the Heaviside operator method or Laplace transform 
method. Recently, another equivalent treatment has been 
given in England by Temple and Jahn employing the method 
of characteristics. In reference 1 a few curves are given for 
the aerodynamic coefficients but no numerical values are 
tabulated. Reference 2 contains no numerical results. 
Temple and Jahn recognize the lack of numerical results and 
supply some initial calculations for the functions necessary 
for flutter calculations. 

A paper has recently appeared by Schwarz (reference d) 
devoted to computing and tabulating the key mathematical 
functions that arise in the theory. The present paper makes 
use of reference 3 to supply more extensive numerical tables 
for application of the theory. The formulas of the theory 
are recast in more familiar form for application to the flutter 
problem and a series of calculations on bending-torsion flutter 
are carried out and discussed. The performance of s imil ar 
calculations for wing-aileron flutter is indicated. Brief dis- 
cussions also are. given of the static instabilities, divergence 
and aileron reversal, and of a one-degree-of-freedom tor- 
sional oscillatory instability. 

For completeness, a connected account of the Possio 
theory is presented since the original presentation in Italian 
is quite terse and also since it is believed that this treatment 
is the simplest and most suitable for general extensions. The 
extension of its application to include the aileron is given. 

AIR FORCES AND MOMENTS ON AN OSCILLATING 
AIRFOIL MOYING AT SUPERSONIC SPEED IN 
TWO-DIMENSIONAL FLOW 

DIFFERENTIAL EQUATION FOR THE VELOCITY POTENTIAL 

The differential equation satisfied by the velocity potential 
in fixed coordinates in the case of infini tesimal disturbances 
is the wave equation 


? ^=V 2 <P 


( 1 ) 


where c is the velocity of sound in the undisturbed medium. 
^For the adiabatic equation of state c 2 =^=y j 
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Figure 1.— Mach angle n. The disturbance at point (£, ij) moving forward with supersonic 
velocity v influences the angular region having half vertex angle g=sln-i?. 

V 


Referred to a system of rectangular coordinates moving 
forward at a constant supersonic speed v in the negative 
sr-direction, the wave equation satisfied by the velocity 
potential in two-dimensional flow becomes 


1 | b 2 <f> | r/tiV , ~j b 2 <j> dV n 

c 2 bt 2 'c? bx df _\c / _ da? by 2 


( 2 ) 


It is proposed to treat the effect of a slightly cambered 
thin airfoil moving forward at a supersonic speed v at small 
(zero) angle of attack as that of a distribution of small dis- 
turbances placed along the r-axis and hence to utilize equa- 
tion (2). The velocity components in the x- and y-directions 
relative to the moving airfoil are, respectively, 


and 




which may be considered the additional components to the 
main stream due to the disturbance created by the presence 
of the airfoil. Relative to coordinates fixed in space, the 
velocity components are v-{-v x and v v . 


z-axis. (See fig. 1.) Upstream from this angular region the 
value of <j> o is zero. It follows also that disturbances hi 
the wake need not be considered and the solution to the 
boundary problem may be attempted by a distribution of 
potentials of the type taken along the projection of the 
airfoil on the z-axis. 

A disturbance at (£, n) created at time T is first fell at a 
point (x, y ) after a certain time ti has elapsed. The point 
(x, y) penetrates the wave front of the disturbed region and 
because it is moving at a speed greater than that of the wave 
front it emerges from the disturbed region at a later time r 4 . 
Thus, the duration of this initial disturbance at (z, y) is 
7-2—7-!. (See fig. 2.) The transition at (z, y) from a region 
of quiescence to a region of disturbance and vice versa is 
associated with the vanishing of the denominator hi equation 
(3) . The values of r, and r 2 for a disturbance created on the 
axis 17=0 are thus given by 

M(x- £) T V(*— SY~y\M 2 — 1 ) , 

Tw — c{M*- 1 ) — : — . (4) 


where the minus sign is associated with n and the plus sign 

V 

with ^ and where M=— It may also bo observed that a 

c 


negative quantity under the radical sign in equation (3) is 
to be interpreted as associated with an undisturbed region 
(that is, with 0=0). 


POTENTIAL FOR A DISTRIBUTION OF SOURCES 

The total effect at any point (x, y ) is the sum of the effects 
of disturbances originating between the leading edge £=£> 
and the intersection of the Mach line through (x, y) with the 
£-axis 

i=ii=x—y-y)M 2 —l 


(since only disturbances created forward of the Mach angle 
region can affect (z, y); see fig. 3). 

The total potential at (x, y) at any time t is thus given by 


EFFECT OF A SOURCE 

Equation (2) is linear and solutions are therefore additive. 
An important particular solution of equation (2) having the 
property of a source pulse is 


- 4 (£> v, T ) ... 

' Vc*(f- TY-[x-ti-v{t- T)f— (y- v ) 2 


(3) 


This solution may be considered to give the effect at a point 
(x, y) at time t of a disturbance of magnitude A originating 
at a point (£, 7/) at an earlier time T. The potential <j> Q is 
thus a retarded potential and the elapsod time at (x, y) since 
the creation of the disturbance is r=t—T. 

Unlike the situation for a subsonic flow, for a supersonic 
flow the effect of the disturbance is propagated only down- 
stream; that is, the point being influenced ( x , y) is always 
considered to be aft of the point of disturbance (|, tj). 
Equation (3) is thus valid in the angular region with vertex 
at (|, ij) and bounded by two straight lines making the Mach 
c 1 

angles ±/i=±sin _1 -=±sin -1 jg with respect to the 


«*■ o-rr * * 

(5) 



Figure 2.— Influence of Impulse created at point ({, 0) at lime f- T on a point (z, t0 fixed 
relative to (£, 01 and moving with supersonic velocity t>. (Observe that the disturbance 
Influences the point (r, g) only during the time Interval n-n.) 
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BOUNDARY CONDITION AND STRENGTH OP DISTRIBUTION 

The function A(£, 0, t— t) giving the magnitude of the 
source distribution is now to be determined by the usual 
boundary condition of tangential flow along the airfoil. If 
the ordinate of any point of the mean line defining the airfoil 
is given as y=y m {x, t), the boundary condition may be 
written 

< 6 > 


where w(x, t) thus represents the vertical velocity induced 
by the source distribution in order to realize tangential flow 
at the airfoil boundary. (In the nonstationary case as in the 
stationary case (corresponding to the Ackeret treatment), 
the two surfaces of the airfoil may be considered as acting 
independently of each other. For the purpose of obtaining 
the oscillating forces in the linear treatment it is sufficient, 
however, to consider separately the upper and lower sides 
of only the mean line.) 


The evaluation of as y approaches zero may be readily 


obtained by use of the 
2t=(tj — Ti)COS fl + T 2 +Ti. 
yields 


variable 6 instead of r where 
This substitution in equation (5) 



Tj — Tl 

2 2 


cos d 



By differentiation with regard to y and with the aid of an 
integration by parts 



Figure 4-. — Sketch Illustrating the three degrees of freedom k, a, and 0 of the oscillating 

airfoil. 


Since %i=x-y-\IM 2 — 1, there results in the limit as y ap- 
proaches zero on the positive side the important relation 


or, briefly, 


A(x,t)=-^w{x,t) 


( 7 ) 


For y approaching 0 on the negative side an equal and 
opposite result is obtained and hence the distribution of 
singularities to be utilized to replace the airfoil is of the 
source-sink type. Thus <j> is to be understood in the sub- 
sequent analysis to be prefixed by a ± sign, +Jbr the upper 
side and — for the lower side. 

The total potential for y = 0 may now be expressed by 
means of equations (5) and (7) as 


l) l 



-,”(*» f ~ T > ---&•# 
V(T— Tj)(t 2 — t) 


( 8 ) 


by c* T 




— t " — - jJL r f r ^ sin 2 5 de d$ 
■y t lv 2 —c t c^jAl 2 —lJa Jo dt 


n 



Fiona e 3.— Sketch, showing that only disturbances created forward of the Maeh angle region 
with vertex at (x can affect (i, »). 


where, from equation (4) with y=0, 


and 


_x—% 1 

Tl_ c ili-t-1 

x — £ 1 

Ti ~ c M-l 


APPLICATION TO OSCILLATING AIRFOIL 

The general result given by equation (8) may now be 
applied for definiteness to the case of an airfoil performing 
small sinusoidal oscillations in several degrees of freedom. 
Let the wing undergo the following motions: a motion due to 
displacement h (velocity A) in a vertical direction; a tor- 
sional motion consisting of a turning about x=x 0 with in- 
stantaneous angle of attack a; a rotation of an aileron about 
its hinge at z—Xi with instantaneous aileron angle /5 meas- 
ured with respect to a. (See fig. 4.) 

In accordance with equation (6) the vertical velocity 
at any point a: of the airfoil situated at 0^x^2b (of chord 
26 and leading edge at 3=0) is easily recognized to be 

w(x, t)= — [A-j-tta-Kz— £o)a+®0-{-(*~ x i)$\ (9) 

where the /3-terms are to be interpreted as zero for x<Cxi 
(and where the minus sign is introduced because the vertical 
velocity w is positive upwards whereas the terms within the 
brackets are positive downwards) . 
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It is convenient in treating sinusoidal motion to utilize 
the complex notation 

h=h^e ial \ 

a=a 0 e i “ t S (10) 

where h „, a 0f and /9 0 are complex amplitudes and hence include 
phase angles. 

Since the further analysis is concerned only with ex- 
ponential time variations of the type given in equation (10), 
the function w(£, t—r) occurring in equation (8) is of the 
form w(£)e ta(,-T) , which may also be written for convenience 
as w(%,t)e~ UoT . The potential <j> given by equation (8) 
may now be written as 

4>(x,t)=^- p*=^= j o (11) 

where 

Tfy S l f T! g- f " r J 


Thus, /(£, a;) is a function of the variable r— £ and of two 
parameters M and w or, alternatively, M and k. 

It is desirable to express the potential 4> as the sum of tho 
separate effects due to position and motion of the airfoil as- 
sociated with the individual terms in equation (13). Thus 

t)=4i a -\-4>h J r4>'a- 1 r4>»-\- < t>b (14) 

where 

vfe “X " I% m 

w= yJ=-i \f« m ' x)di 

462 C x 

^ i= VM 2 -“ “J, (f-ab)Af. *)# 

Oh 

FORCES AND MOMENTS 


The integration with regard to r may be readily performed by 
substitution of the variable d where 2r=(r s — ti) cos 0 +T 2 +T!. 
Then 


a;)=— g-fu<rri-n)/2 f g-b 00S4 (ti-t0/2 


With ri and t 2 replaced by their values as given for equation 
(8) and with the aid of the Bessel function relation 


I (V* coa *cf0=:«/o(X) 


it is recognized that 


/(£,»:) = e 


Af 

c M>-1 T 
«0 


( &Z. I M \ 

V c M*~ l) 


( 12 ) 


The basic pressure formula in the theory of small dis- 
turbances is 


P= 



which in the present case of the moving airfoil may be 
expressed as 


P=-p 



b<f>\ 

v bx) 


where p is the density in the undisturbed medium. The local 
pressure difference on the airfoil surface between the upper 
and lower surfaces at any point x (nondimensional) is 


p'=-2p 


c )4>\ 

\bt~'2b dr/ 


(15) 


Throughout the subsequent analysis it is convenient to 
employ the variables x and £ in a new sense to mean non- 
dimensional quantities obtained by dividing the old variables 
by the chord 2b. The retaining of the symbols x and £ for 
the nondimensional variables should lead to no confusion. 
The potential <t> of equation (11) is then 

4>0M) = V A/^-l Jo (£ — ; Zi)/3]/(£,z)d£ 

(13) 

where with the introduction of theimportantfroquencyparam- 
eters 

k=^ 

v 


The total force (positive downward) on the airfoil is 
P= 2 b£ p'dx 

— ^- 2 pyj^ ^dx— 4 pb j <t> dx ( 16 ) 

The moment (positive clockwise ; fig. 4) on the entiro 
airfoil about any point Xq is 

M 0 = 46 2 J' (x—Zo)p'dx 

=— 4 pbvj^ ^ (x— Xz)dx— 8 p& 2 J o 4>{x—Xo)dx ( 17 ) 


Cxi 


2kM 2 


the function /(£, x) becomes 


/(£, x)^e-^-vj 0 [^i)] 


(120 


Similarly, the moment (positive clockwise; fig. 4) on the 
aileron about the hinge point x t is 



=— 4 pbvj ^ (x—xi) dx—8pb 2 j 4(x—X!)dx (18) 
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In the further reduction of equations (16) to (18), with 
the potential <f> replaced by its separated form given in 
equation (14), the following sets of integral evaluations are 
required: 

j; dx=- r ^_- a [r i (M,k)—x 0 r 1 ( M,k )] 

C 1 2b 

J. i*=-j w = l mr,W,k) 

( *;*— a k)-v,(U, t)J 

ri 2b 

J (ftp dx 7 — 2 ky £1) 

£ 4> t i dx = i $ | ; j( m > k > 

X' ^ 1 *- vfe 4 B ^ «] 

J" ^ x dx= lyt—i C<3 x i)l 

£ ^ x dx= ^L- |S [| s s (M, A, x.) +x,f 2 (M, A, *)] 
ij <f>a x dx— ~[^£i | 0® 2 2*(^i 1) 

£ 4>aX dx= ^|==== a jji ffj(M, A) — | Xo g a (M, A)J 

p <&rx dx= x v & £§ s s(^i *> x i)+ *i) J 

J <^x dx=-^==^ $ s a (M, 1’, xj xMM, k, X,)] 

£ (2-^)*=^= 0a Pi(M, A", XO 

£ ^ (X— x 1 )tfx= -^=== ^ d |^| p»(M, A, Xj)— XoP^M, A, Xj)J 
X l)x (x— x I )dx=-jJ=^ 00 s 1 (M,A,x 1 ) 

X ^ * 5 S ’ (M ' *■ *> 


X i 26 1 

( j >a (x—x 1 )dx= -jjj=j va ^ Ps(M, A, xj 
£ 4 >i(x—x 1 )dx= ^==y d ^ Ps(M, Ar, x x )— ^ x 0 p a (M, A, xoj 

X i 26 1 

<f> 0 (x-x 1 )dx= , Mi _ £ ®/3 2 S *( M > *» *i) 

f 1 46 2 • 1 '■ 

<t>i,(x—x 1 )dx=-jj = == = /3 g s 3 (M, A, x x ) 

The functions defined by the foregoing integral evaluations 
are further discussed in the following section ; first, however, 
the force and moments (equations (16) to (18)) are given in 
their final forms as 

P= tJ=== pKto-f- X— 26x 0 a)r 1 +26(2m+?i— 26x 0 a)r 2 + 

1 L 

46 2 d |+r 2 /3f 1 +46^f 2 +46^ |] (16') 

M«=~ £g(ro!+A— 26xb«)gi + 26(2pd+A — 26xo«)^+ 

46 2 d |+» 2 j3(s 1 +x 1 f l )+46^ (| 2 +x,f 2 )+ 
46 2 £}(|+x 1 |)]-26x 0 P (170 

p (m+A— 26xod) 2^+26 (2rd4-A— 26xo«) 

46 2 d|+0 a /3s 1 +46^|+46 2 /3|] (18') 

REDUCTION AND EVALUATION OF FOREGOING INTEGRALS 

It is convenient to introduce the substitution u=x— £ and 
to express the function /(£, x) (equation (12')) as 

I{Z,x)=I{u)=e-“*J,(j £ u) (19) 


The various functions defined by the foregoing sets of 
integrals may now be expressed as follows: 


n {M, k ) = £ / (u) du 
r,(M,» = XX I (u) du dx 
r.(M,«-2XX (x— u) I (u) du dx 
g, (M, k) = £ uI (u) du 

fc (j/,e-aXX x7 (u) du dx 

fc^tt^xx x (x— u) J(u) du dx 
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k, x x ) 
Pi(M, k, x x ) 
p s (M, k, x : ) 
ti(M, k, Xi)-- 
U{M, k, 2!) = 
ti(M, k, x,)-- 
S\{M, k, x{)-- 
8 2 (M, k, x x ) = 
Si{M, k, x x )= 


= f (u— x x )I(u)du 

JX i 

= 2 J J* (x—x x )I{u)du dx 
= 6 f f (x—Xi)(x—u)I(u)du dx 

Jx i JO 
ri-xi 

= I I(u)du 

r i-ii rx 

■J, J. 

-rr (x— u)I{u)du dx 

ri-ii 

= I ul(u)du 

pi-xi fx 

-2 J J xl(u)du dx 

X i -a rx 

I x(x—u)I(u)du dx 


Borbely (reference 2) has shown by means of reduction 
formulas that the six r- and ^-functions may be obtained 
from a single integral. In a similar manner it may be indi- 
cated how the foregoing 16 functions may be obtained from 
the evaluation of the same integral. The reduction is 
accomplished in two stages. First, consider integrals of the 
following type: 


h=(l— x x )K 
fa=(l— h) 

t%= (1 — Xi) 8 (Ao — 2&i -j- hi) 

®1 = (1 — 2l) S ^I 

8 2 ={\—X l )\h Q —h i ) 

Sg— ( 1 — x x )*(2ho — 3Ai + kg) 

The final stage in the reduction of these functions is to 
utilize the following recursion formula (reference 2) obtained 
by integration by parts: 

3) «[*+ ( 1 - X) IJ-'V.Qj) - + 

i(l-2X)/x_ 1 (M,5) + 

(1-X) 2 |/x_ 2 (M,5) (21) 

where X^ 1 and/ with a negative subscript is to bo interpreted 

/ M 2 — 1 _ \ 

as zero. (Observe that ■ - ^p u>=2k.\ 

The function /x (M, 5) may clearly refer also to the fore- 
going g- and A-functions, if 5 is replaced by tho appropriate 
parameter; namely, Sx x for g\ and 5(1 — x x ) for h\. (See 
equations (20).) The recursion relation (equation (21)) 
thus reduces the various functions to the single function 

/„ (M,«)=iJV“ J„ Q) du ( 22 ) 


/x=/x(M, co)=J^ I(u)v>du j 

0x=/x(M, wr 1 ) = -4pi f I(u)uHu } (20) 

*i Jo l 

_ j ri-x, | 

h\=j\[M, w(l— x 1 )] = ( 1 — ari )x+i J 0 I{u)vkdu J 

By integration by parts it can be readily verified that the 
following relations hold: 

r i—fo 

ra=fo~fi 
rt=fo— 2/i+/s 


which is therefore the only integral needed in the evaluation 
of the forces and moments. 

The important integral in equation (22) has been recently 
made the subject of a mathematical investigation by Schwarz 
(reference 3). Schwarz gives tables of the values of its real 
and imaginary parts to eight decimal places for OS «S5 
and for IS MS 10 for conveniently small intervals. For 
values of 5j>5 not given in Schwarz’ tables, the function / 0 
may be evaluated by means of the following series develop- 
mept (reference 2): 

f 0 (M, «) ") 2 "n! ( 2 n+l ) [Jn 

( 23 ) 


2i=/i 

2s=/o — ft 
g» — 2/o — 3/! +/ j 
p l =q l —x l r 1 +x 1 2 {g a —g 1 ) 

Pi=q 2 — 2x^3+ x? (g 0 — 2g t + g t ) 
Pi=ii— Sx^+xSigo- Zgi+Zgi—gi) 


Table I gives values of the functions / 0 (M, w) based on 
the tables of Schwarz and on equation (23) for selected 

values of the Mach number M=~> y> ~> 2> |> -^» and 5 


and for various appropriate values of w 



Later use 


is made of the values given in table I for obtaining tables 
for flutter calculations. 
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EQUATIONS OF MOTION AND DETERMINANTAL EQUATION 
FOR FLUTTER CONDITION 

The equations of motion and the border-line condition of 
unstable equilibrium yielding the flutter speed and frequency 
may be obtained exactly as in the incompressible case treated, 
for example, in reference 4. The two-dimensional treatment 
(infinite aspect ratio) is retained herein. Modifications due 
to assumed vibration modes of the finite wing may of course 
be introduced as in current practice (for example, reference 5). 
The modification of the forces and moments due to the three- 
dimensional nature of the flow is a more difficult problem 
which remains to be studied. 

The equilibrium of the vertical forces, of the moments 
about the torsional axis x=x a , and of the moments on the 
aileron about its hinge x—x 1 yields the three equations 


AM-4- aiSa-H fiSf} -f- A C ,,= P 

6cl a ~}r 3[lis-f-26(iCi — r-o)S^] -j~ hS a -{-aC a =Ala 

j3li}-\-'a[Ip-\-2b{xi — xo)S e ]+hSe+pCe=M fl l 


(24) 


where the various parameters are defined in the list of 
notation. (See appendix.) 

In order to define the borderline condition of unstable 
equilibrium separating damped and ■undamped oscillations, 
the variables A, a, and /3 are used in the sinusoidal exponential 
form given in equation (10). For the desired condition, it is 
necessary that the equations (24) have a (nontrivial) solu- 


tion for the complex amplitudes A 0 , a 0 , and /So, or that the 
following determinants! equation hold: 



A ca 

Acp 



A a ji 

A aa 

Aafi 

= 0 

(25) 

Am 

A. ba 

At# 




where the complex elements of the determinant in separated 
form are 

Ac = — jLix a -f- Lz + i L t 

Ac$= — nxp-{- 

Aan— —itXa+MiAiMj 

A aa =Q a X — pir a 2 -\- j 1/j -f- 

A a f=—p.[r fi 2 +2( i x 1 —x ii )Xff]+M i +iAI t 

■d 4 j,= — jur/j+iVi+iA^ 


A ia = — /t[r 3 *+2(®i— xbJzrf+iVi+tiV'* 
Atf— QpX — A/j-f-iA/j 


and where the X’s, M’s, and A T ’s are defined by the force 
and moment equations (16'), (17'), and (18') expressed in 
the following forms: 


Hence, 


M a -- 

Mp ■■ 


t— N ' — 

-4 p&VAV"' j~ (Mj+iM*) + ao(Mj+iA/*) + 0o(M*+iM s ) 
-4pb i c i k t e , * t [(^) (JVi + iN s )+ a 0 (N 3 +iNJ + A(M+W,)] J 

;(- 2 ^+r0 


Li+iLt= ^m=i x 

X 8 +ii 4 = = [-2r*+f r s -£ (^-2r,+j r,)- ^ 


-2*o(-2r s +£n)] 


Mi+iMi = {"1 2a+ f 

M s +iM a = 2_ j { — 3 


M+iM 2 =^p =T [- 223+1 a-2*(-2i*+J*)] 

(~ 2gj+-j- 2i^~ — 2r s -(— ^-r 3 — ^ ^ 2^+^^— 2g 2 +-j- 2r 0 

f *- j (~2s*+J *i)+2(zi-Zo) [-2^+| (~2^+| ti) 

XT I -A7 1 I 2t _ \ 


/ I— 

Ni+m== jAP-i ( _2?2+ r Pi ) 

Ni+iNi= jm=! [ _ i Pi ~i ir 2pi+ T P')~ 2x ° (- 2pi+ T pi )~] 

1 T 4 . 2i i( 0 . 2i 

= 3 8s+ J 8t ~k V _2s3+ ¥ 


Ns+iNs 


2 % \ 
~T U 
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The detenninantal equation (25) with the foregoing com- 
plex elements is equivalent to two real simultaneous equations 
and hence may be solved for two unknowns. In a given case 
the usual unknowns are the flutter speed v and the flutter 
frequency u or, more conveniently, the related nondimen- 
sional parameters X and 1 jk. The parameter X appears 
linearly and only in the major diagonal elements (with bars) , 
while the parameter 1 jk appears transcendentally in every 
element of the determinant. Hence an obvious procedure, 
though not the simplest for obtaining the simultaneous solu- 
tions of the two equations, is to fix values of 1/k, to solve for 
the roots of the two polynomials in X, to plot graphically 
these roots against 1 jk, and to note the points of intersection. 

In a systematic numerical study of flutter any two param- 
eters may be utilized as unknowns instead of X and 1 jk, a 
procedure which is often more convenient. A discussion of 
such procedure and the use of a method of elimination for 
simplifying the calculations is given in the appendix of 
reference 6. 

The application to the two-degree-of-freedom subcase of 
bending-torsion flutter is treated more fully in the following 
section. 

APPLICATION TO BENDING-TORSION FLUTTER 

The determinantal equation in the two degrees of freedom 
h and a is 




X ca 

■A-aa 


= 0 


or 


QhX — p-\-Li-\-iLi 
— itXa + Mi + \M% 


U a X— pr a *+M s + iMi 


=0 (27) 


The two equations in X obtairied by equating the real and 
imaginary parts separately to zero are 


Q.j£l a X 2 + [fia(£l — A — jJ.r a a )]X- j- C R =0 

(Oaij-f-QiAf^AT-)- Cz= 0 

where 

C R =p[x tt (Mi+JU) — (M 3 —pr a 2 ) — InrJ— pxj] +D t 
C I =p[x a (M i +Li ) — Mi — LirJ] -\-Di 
and where 


(27') 


D R =I^Mi — L^Mi — LiMi-\- LiMi 


Dj — LiMi — LiMi -f- T13M1 — ijA/j 


For convenience in numerical tabulation, it is desirable to 
introduce primed quantities, independent of the parameter 


x 0 , defined by the following relations: 

— 2xJLi 

Li = Li — 2ro L<z 

M 1 =M 1 '—2x 0 L 1 

2xoLi 

M 9 =M 3 '-2xoKM 1 '+L a t )-2x Q L 1 )\ 

il/,= Af 4 '-2ro [(MZ+L^-ixM ) 

In table II convenient expressions for the quantities L u 
L 2) Lt , Li, Mi', M 2 ', Ms, and Mi are given and tabulated 
together with the combinations M/+L3' and Mi'+L/. 
Clearly these quantities depend on the function / 0 given in 
table I and hence the tabulation is made for the same values 
of M and 1 fk (or «). In addition, table II contains values 
for the quantities D R and D t which, in fact, are independent 
of Xo and may be expressed as 

D r = LiM z ' - L t 'Mi ' - Li Mi' + L^ M t ' 

D r =LiMi' - L t ' Mi' + L 2 M/-L t 'Mi' 

The numerical application in the case of bending-torsion 
flutter has been performed for various selected examples. 
In most of the calculations the numerical procedure was to 
fix values of 1 jk, eliminate X, and solve for the parameter 
x a . Interpolation was also used to obtain additional points 
in order to improve the failing of some of the curves. Values 
of 1 jk less than 1 did not yield any flutter points in this 
procedure. Results are shown plotted in a number of 
figures (figs. 5 to 20); however, before these figures aro 
discussed, it is desirable to explain the significance of the 
parameters and the numerical values assigned to them. 

The parameter p may be considered to signify the wing 
density and three selected values 3.927, 7.854, and 15.708 
in the order of increasing wing density have bceu mainly 
used in the calculations. (These values correspond to 

values of j-=5, 10, and 20 in the notation of reference 4.) 

Alternatively, an increase in p may be interpreted as an 
increase in altitude for a fixed wing density. The parameter 
P may be expected to range up to high values for actual 
supersonic wings at high altitude. Only a few calculations, 
however, have been made for high values of p (p= 78.54, 

^=100; see fig. 18). 

The parameter is the ratio of the wing bending 

frequency to the wing torsional frequency and may be 
expected normally to be less than unity. The three values 
0, 0.707, and 1 have been largely used in the calculations 
although other values up to 2 have also been studied. 

The parameter x 0 represents the position of the elastic axis 
measured from the leading edge and the three values 0.4, 
0.5, and 0.6 represent, respectively, positions at 40, 50, and 
60 percent chord. (These values correspond to values of 
a=— 0.2, 0, and 0.2 in the notation of reference 4.) 


’ ' 

( 28 ) 
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The parameter x a represents the distance of the center 
of gravity from the elastic axis. For example, x a =0.2 
represents a position of the center of gravity 10 percent of 
the chord behind the elastic axis. In many of the calcula- 
tions x a has been regarded as variable. 

The parameter rj represen ts the radius of gyration of the 
wing about the elastic axis and has been kept fixed at the 
value r a s =0.25. 

The ordinate in figures 5 to 20 is the nondimensional 
flutter coefficient v/bo} a where bu a is a convenient reference 
speed. This coefficient is also a function of the Mach 

number .lf=~ and several values of M have been employed 

in the calculations. 

In a plot of the flutter coefficient v/bu a against M, straight 
lines drawn from the origin at angle 5 and intersecting the 
curves may be given an interesting interpretation (fig. 17). 

The slope of the line is given by or cot & = ~~" 

Thus, cot 3 is directly proportional to the product of the 
chord and the torsional frequency divided by the velocity 
of sound. The question of whether at a given value of M 
the value of boi a which will just prevent flutter is also suffi- 
cient to prevent flutter at neighboring higher values of M is 
answered by the simple criterion of whether cot 8 increases 
or decreases. In figure 17 two typical flutter curves are 
shown. In curve B the value of bu a just necessary to pre- 
vent flutter at a speed corresponding to the value of M 
at P 2 is insufficient to prevent flutter at any higher value 
of M for which the flutter curve is below the straight line 
OP*. For the type of curve A a maximum value of 3 occurs 
at the “design critical points” Pi. The value of bu a just 
necessary to prevent flutter at a speed corresponding to 
the value of M at Pi is also sufficient to prevent flutter 
at all higher speeds. 

The. following salient facts may be extracted by inspection 
of the figures. Flutter exists or is possible for various 
ranges of the parameters but, in general, compared with 
subsonic cases the ranges of the parameters yielding flutter 
are more restricted. 

The chordwise position of the aerodynamic center, the 
center of the oscillating pressure, is an important factor in 
the consideration of flutter. In the static case the mid- 
chord is the aerodynamic center for AT>>1. For subsonic 
speeds, M<<jl, the linearized theory indicates the quarter- 
chord position as the aerodynamic center. It should be 
expected that in the transonic region near M= 1 the aero- 
dynamic center may shift considerably. From this point 
of view alone conclusions drawn from the simple theory 
for the range near M= 1 may require large modifications. 


The nature of the modifications may be roughly inferred by 
further experimental and theoretical study of the behavior 
of center-of-pressure locations. 

For low values of the ratio of bending frequency to tor- 
sional frequency — »0 the position of the center of gravity 

Cx ) a 

relative to the aerodynamic center is important. For center- 
of-gravity positions forward of the midchord no flutter 
exists, whereas for positions behind the midchord there is a 
sharp decrease in the flutter coefficient from infinity; the 
position of the elastic axis influences the value of the flutter- 
coefficient in this region, forward positions being more fa- 
vorable (figs. 5 (a) to 16 (a)). 

For values of the position of the center of gravity 

Cl3 a 

relative to the elastic axis becomes of more importance. 
For center-of-gravity positions forward of the elastic axis 
no flutter exists, whereas for positions behind the elastic 
axis flutter does occur, and a relative minimum coefficient 
appears for center-of-gravity positions only slightly (a few 
percent of the chord) behind the elastic axis. 

The intermediate case, for which —=0.707, shows a blend- 

C0<x 

ing of the effects in which the center-of-gravity position 
relative both to the aerodynamic center and to the elastic 
axis is significant. 

In figures 12 and 14 there are shown, for reference, some 
numerical values of co/w a , the ratio of the flutter frequency 
to the torsional frequency. 

The effect of the wing density parameter n is rather com- 
plicated but, in general, an increase in ^ yields a correspond- 
ing increase in the flutter coefficient. For low values of 
w*/&! a and for high wing densities this increase is expected to 
be proportional to -fy.. In the resonance-like region near 

— =1 and for small values of x a the flutter coefficient is rela- 

C0 a 

tively unaffected by the value of fi, and in this region the 
structural damping may be expected to be particularly 
effective in increasing the flutter coefficient. 

For values of the Mach number near unity (for example 

AI=^> a value for which the validity of the theory is in 

question), the flutter calculations become difficult to plot 
because of the appearance of other branches. In some cases 
(for instance, £6=0.6) the flutter instability appears limited 
to a definite range of flutter speed coefficients. Calculations 
to include damping were performed to verify the existence of 
the range. (The appearance of these other branches seems 
to involve values of 1/k for which the quantity ili* is negative. 
The condition of negative il /« is significant for the one-degree- 
of-freedom instability discussed in the next section.) 
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(a)^-O. • tt>) ;r -0.707. Co) ^-l.o. 

Wh Ua (da 


Figure 9.— The flutter coefficient against center-of-gravlty location for several posltfons of elastic axis and for three values of the frequency ratio. A/— tt: ^-7^64. 



(a) 77-0. (b)— -0.707. (c)— — 1.0. 


Figure 10.— The flutter coefficient against center-of-gravity location for several positions of elastic axis and for three values of the frequency ratio. Af — y;>i-lfr.708. 
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FlOVSE IX. — The flutter coeffldent against center-of-gravity location for several positions of elastic axis and for three values of the frequency ratio. M-1; jt-3.927. 


(b) — -0.707. 
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(o)— -1.0. 
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measured from leading edge, 
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Center— of -gravity location measured from elastic axis, percent chord 


(b) ^-0.707. 
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(o) ”*“1.0. 
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Figure 12.— The flutter coefficient against center-of-graviiy location for several positions of clastic ails and for three values of the frequency ratio. M—2; p- 7.854. 
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Center-of-gravity location measured from elastic axis , percent chord 


(b) — -0.707. 
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Figure 18.— The flutter coefficient against center-of-gravity location for several positions of elastic axis and for throe values of the frequency ratio. Af-2; ^-15.708. 



Cents/ — of- gravity location 
measured from leading edge, 
percent chord 


Center -of -gravity location measured from e/astio axis, percent chord 


(a) o. 


(b) — -0.707. 


Figure U.— The flutter coefficient against center-of-gravity location for several positions of elastic axis and for three values of the frequency ratio. M-S\ g- 3,927. 
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(a)— -0. 


0»E-O.707. 


Ud 


FrGVBE 16.— The flutter coefficient against center-of-gravity location for several positions of elastic axis and for three values of the frequency ratio. M-5; *-7.854. 



(a)— — 0. (b)^— 0.707. , (c)^-lA. 

Uc COm Wa 


Figure 16.— The flutter coefficient against center-of-gravity location for several positions of elastic ails and for three values of the frequency ratfo. Af— 5; *—15.708. 
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FIQcbk 17.— The flutter coefficient against Macb number for two locations of the center of 

gravity. Other parameters are ——0.707; a-0: u-7.854. 

<*> 

A plot of the flutter coefficient against Mach number for 
two values of x a is shown in figure 17. The significance, of 
the straight lines drawn from the origin has already been dis- 
cussed. The type of curve A is representative of tlie effect 
of forward location of the center of gravity and the type of 
curve B is representative of rearward locations of the center 
of gravity. Figure 18 gives a plot of the flutter coefficient 
against M for various values of the wing density parameter p 
and for a rearward location of the center of gravity. The 
subsonic values for M= 0 and M=0.7 shown on these curves 
and on some of the other figures have been either taken from 
reference 7 or calculated in the manner outlined therein. 
The subsonic and supersonic parts of the curves (figs. 17 and 
18) have been arbitrarily joined by a dashed smooth curve in 
the transonic range. In figure 19 there is given a cross plot 
of flutter coefficient against frequency ratio «„/w a , for various 
values of M, and in figure 20 is given a similar cross plot for 
three values of the elastic-axis parameter x 0 . 

An indication of the effect of structural damping in in- 
creasing the flutter speed in a few examples may be obtained 
from the following table, where g a and g h are the torsional 
and flexural damping coefficients, respectively, and where 

M=-y)^=7.854, a=0, and a-„=0.2: 


OhJcOa 


h 

«/«* 

P/&W* 

0 

0 

0 

0.673 

2.438 

0 

.06 

0 

.M3 

2.551 

0 

.10 

0 

.628 

2.669 

.707 

0 

0 

.777 

1.535 

.707 

.00 

0 

.771 

L 653 

.707 

.10 

0 

.766 

1.669 

.707 

0 

.06 

• .788 

L 692 

.707 

0 

.10 

.797 

1.642 

.707 

.06 

.06 

.782 

1. 623 

.707 

.10 

.10 

.784 

1.725 


STATIC CASES— WING.DIVERGENCE AND AILERON REVERSAL 

It is of some interest to examine the expressions for the 
forces and moments in the limit case in which the frequency 
approaches zero. There follow then for the mean-line wing 
section the well-known static-ease results which may of 
course be obtained directly without the use of a limiting 
process, as originally treated by Ackeret. Thus, with the use 
of the following relation easily verified from equations (20), 

Iim/x (m, k) =r4-r 

fc - X > A -+-1 

there are obtained from equations (16') to (18') for the lift 



Fiovri 18.— Tho flutter coefficient against Mach number for several values of p. Otlicr 
parameters are — -0;r«-0.2;a-0. 
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and moments in the static ease, 




Ma 


4 




[(1 — 2a%) ce-J-(l — *) (1+an— 2^b)j8] 




These relations for the mean-line wing section are now 
used to obtain the critical speeds for the static instabilities — 
wing divergence and wing-aileron reversal (for wing of 
infinite span). At the wing divergence speed the effective 
torsional stiffness of the wing vanishes, hence the total 
moment about the elastic axis is zero. The sum of the struc- 
tural restoring moment and the aerodyn ami c twisting 
moment is 


aC a 4 


4 p6V 


( 1 - 2 *,) 



Fksuie 19— Tbc flutter coefficient against frequency ratio for several values of M. Other 
parameters are c=0; r«=»0.2; 7.8J4. 



parameters are jtf-ylj.-O.fcjt-ZSM. 

which when equated to zero yields the divergence speed 

Thus, the divergence speed is real only for positions of the 
elastic axis behind the aerodynamic center (midchord, in the 
simple theory) . This formula obviously should not be used 
for values of AT too near unity. 

For comparison it is of interest to note the corresponding 
result for the divergence speed in the subsonic case, where the 
aerodynamic center is (approx.) at the quarter-chord point. 
Thus, 


where M is less than about 0.7. 

The aileron reversal speed is determined by the condition 
that the change in angle of attack due to wing torsion nullifies 
the effect of movement of the aileron so as to yield no change 
in lift (in rolling moment, in the case of finite wing span). 
There are two equations to be satisfied for this condition; 
namely, 

«+(i— *00=0 

(that is, L= 0) and 

*C a + [(1— 2*,)«+(l-aO(l+* 1 — 2*,)fl=0 
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Tlie aileron reversal speed, obtained by elimination of a 
and /3, is 

v R =bo> a (M 2 — ly^^prj-jzr 

yxi 

For hinge positions aft of the midchord, the factor l/ya^ in 
this expression varies from 1.4 to 1.0. The aileron reversal 
speed is thus relatively unaffected by the position of the 
hinge. In general v R may be expected to be lower than v D . 


ONE-DEGREE-OF-FREEDOM OSCILLATORY INSTABILITY 


As was pointed out by Possio, the theory indicates the 
existence of a torsional instability which may arise for a wing 
having only one degree of freedom. This instability is due 
to the wing being negatively damped in torsion and is asso- 
ciated with the vanishing (and change in sign) of the tor- 
sional damping coefficient M* (equation (26)). 

Certain considerations for the case of slow oscillations 
made by Possio (reference 1) and further discussed by 
Temple and Jahn serve to bring out the main results. Thus 
from equation (20), for slow oscillations, 


and 




1 

x+i 


. 2 kM 2 1 




The condition M,(M, x 0 )— 0 is shown plotted in figure 21, 
where the shaded area is the region in which the instability 
is possible (negative M t ). The maximum ranges for the 
parameters r 0 and M in this region are x a less than 2/3 and 
M less than V2h (and greater than Unity). 

(It may be appropriate to mention that a similar torsional 
instability is theoretically indicated even in the subsonic 
(incompressible) case for positions of the axis of rotation 
between the leading edge and the quarter-chord point. The 
combination of parameters required for tills indicated insta- 
bility, however, is not very likely.) 

The torsional instability may be studied more fully in the 
general case. It is found that the range of instability for 
the parameters x 0 and M remains essentially as in the simple 
case (large 1/k) but more information may be obtained re- 



garding the critical speed and frequency. The moment 
equation is equivalent to A aa =0, or to the two equations 

nrj +M S (M, Zo) = 0 

Mi(M, Xo) + g„!U=0 

where the structural damping coefficient hi torsion g„ has 
been introduced as in reference 6. The critical speed and 
frequency may be studied as functions of the parameters x 0 , 
M, g a and the product combination ixr a *. Results of a few 
selected calculations are shown plotted in figure 22. Since 
instabilities are indicated for the range of near-sonic values 
(1<M<1.58), it would seem that a more comprehensive 
investigation of this problem is very desirable. 

It may be remarked that a similar analysis for pure 
bending exhibits no instability while the case of the aileron 
alone does exhibit a range where such instability may 
occur. _ This range for an aileron hinged at its leading edge 
is 1<M^V2. 


Langley Memorial Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., May £9, 1946. 
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(a) Flutter coefficient against aite-of-rotatlon posi- 
tion lor several values of -V (a-16.708). Note that 
the range of i» narrows with increase in AT and 
disappears at Af-1.68 and ra= 0.33. 


(b) Flatter coefficient against torsional damping co- 
efficient Tor two values of a 18.708^ - 

Negative damping values are shown dashed and 
have no physical existence. 


(e) Flutter coefficient against wing density param- 
eter it for several values ofii^Af— The 
stralght-llna curve shown corresponds to A/i— 0 
fo=0.327). 


Figure 22.— Curves for one-degree-of-freedom torsional Instability. 


APPENDIX 


SYMBOLS 


$ 

t 

T 

r=t — T 


P 

7 

V 

c 

M 

X 

y 

a-o 




disturbance velocity potential 

time at which disturbance influence is felt 

time at which disturbance is created 

pressure 

pressure difference 
density 

adiabatic index (for air, 1.4) 
velocity of main stream (supersonic) 
velocity of sound in undisturbed medium 
Mach number (v/c) 

coordinate measured in direction of main stream 
ordinate 

abscissa of axis of rotation of wing section (elastic 
axis) 

abscissa of aileron hinge 

abscissa and ordinate of point of disturbance 


b one-half chord 


After equation (12) the quantities x, y, x 0 , x u and £ are 
employed nondimensionally and are referred to the chord 2b 
as reference length. 


w(x, t ) 

h 

a 

P 


k 


Hi, *) 

JnW 


vertical velocity at position x on chord and at 
time t 

vertical displacement of axis of rotation 
angular displacement about axis of rotation 
angular displacement of aileron; measured with 
respect to or 

angular frequency of oscillation 
reduced frequency (ub/v) 


frequency parameter 



function given in equations (12) and (12') 
Bessel function of order n 
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The following additional symbols, ejnployed in the flutter 
equations, conform to the notation used in references 4 and 
6, in which the half-chord b is the unit reference length. 


M 

S a 

Sg 

la 

h 

a 

c 

x a 


Xg 


r a 


rg 


C„ 

Gg 

Gg 

U a 


03g 




mass of wing per unit span 
static moment of wing-aileron combination per 
unit span referred to elastic axis 
static moment of aileron per unit span referred to 
aileron hinge 

moment of inertia of wing-aileron combination 
about elastic axis per imit span 
moment of inertia of aileron about its hinge per 
unit span 

coordinate of elastic axis measured from mid- 
chord (2xq— 1) 

coordinate of aileron hinge axis measured from 
midchord (2z x — 1) 

location of center of gravity of wing-aileron system 
measured from elastic axis S a /AI b ; location of 
center of gravity in percent total chord measured 

from leading edge is 100 (^b+T^ 

reduced location of center of gravity of aileron re- 
ferred to c ( Sg/M b ) 

radius of gyration of wing-aileron combination re- 


ferred to a 



reduced radius of gyration of aileron referred to c 

(VS) 


torsional stiffness of wing around a per unit span 
torsional stiffness of aileron system around c per 
unit span ___ 

stiffness of wing in deflection 
natural angular frequency of torsional vibrations 


about elastic axis 



(w a =2ir/a, where j a 


is in cycles per second) 

natural angular frequency of torsional vibrations 
of aileron around c 
natural angular frequency of wing in deflection 




/x 1 1/ \ 

wing density parameter (7 - or 4^52 ) (Note that 

in the incompressible case (references 4 and 6) 
,u is replaced by 1 /k.) 

ratio of mass of cylinder of air of diameter equal 
to chord of wing to mass of wing, both taken for 

equal length along span (This ratio may 

be expressed as «=0.24 Off) (j~) where If is 


weight in pounds per foot span, b is in feet, and 
p/po is ratio of air density at altitude to that for 
normal standard air.) 

9a, 9e, Qh structural damping coefficients (see reference 6) 
Li, Li, Lz, Li, Mi, M 2 , Mg, Mi quantities defined in table II 
and by equations (26) and (28) 
v/bu a flutter coefficient; that is, flutter speed divided by 
reference speed bu a 


*•?(?)' 

X— iir a 2 0~) for case of bending torsion 
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TABLE I —VALUES OF / 0 (M, 5) = (/,) s +i(/c)j 


< /4)£ "iJo 3j4 (s) oos “‘ i “ 

J7 170 (h) stnu<i “ 


5 

1 

* 

w* j 

( ft ) i 

03 

1 

ifc 

oa* 

(« i 

5 

1 

k 

( a >« 

( jt)i 

Of 

1 

* 

(A). 

(ft)/ 

*-'{ 



3 



20.00 

12.00 

iaoo 

9.00 

8.00 

7.00 

6.00 
6.00 

4.20 
8.60 

3.20 
2.80 
2.60 

2.20 
2.10 
1.90 
1.68 
1.50 
1.40 
1.80 
1. 10 
LOO 

.90 

.84 

.80 

.70 

.00 

.66 

.62 

.36 

.26 

.16 

.10 

0. 826.12 
.87719 
1. 05263 
1.16959 
1.31579 
1.60376 

I. 75439 
2.10628 
2.50627 
2.92398 
3.28947 
3.76940 
4.21053 
4.78409 
6.01253 
5.54017 
6.29560 
7.01764 
7.51880 
8.09717 
9.56S3S 

10.5263 

II. 6959 
12.6313 
13.1579 
15.0376 
17.6439 
18.7970 
20.2429 
29.2398 
40.4858 
65.7895 

105.263 

IjB 

b 

B ujs.i 

1 

■ . 

—0.14999 
-.21130 
-.21774 
-.23999 
-.24647 
- 24545 
-.27930 
—.30624 
-.29526 
-.30153 
— . 12996 
-.37017 
—.42388 
46577 
-.47752 
-.49523 
-.50275 
-.49720 
—.45896 
-.47679 
-.44031 
-.41597 
-.38765 
-.30881 
-.36551 
-.81977 
-.28073 
-.26126 
-.24735 
17576 
-.12839 
-.07962 
-.04991 

20.00 

iaoo 

6.00 

4.40 
3.30 
2.80 
2.60 

2.40 
2.20 
2.00 
1.96 
L84 
1.66 
1.66 
1.48 
L34 
L24 
1.10 
1.06 

.98 

.94 

.88 

.82 

.78 

.74 

.70 

.66 

.60 

.56 

.52 

.48 

.44 

.40 

.36 

.34 

.32 

.30 

.28 

.14 

.00 

0-27778 
.56556 
1.11111 
1.26263 
1.68350 
1.98413 
2. 13675 
2.31481 
2.52625 
2.77778 
2.83447 
3.01932 
3.34672 
3.66126 
3.76375 
4.14594 
4.48029 
6.05051 
5.24109 
5.06893 
5. 91017 
6.31313 
6.77507 
7.12251 
7.60751 
7.03651 
8.41751 
9.25926 
9.92063 
10. 0838 
1L 5741 
12.6263 
13.8889 
15.4321 
16.3399 
17.3611 
18.5185 
19.8413 
39.0825 
92.5926 

-a 02589 
.02530 
.19004 
.21539 
.22644 
.26381 
.29327 
.33266 
.38184 
.44053 
.45329 
.49337 
.55768 
.69498 
.62531 
.67878 
.71672 
.76856 
.78294 
.81095 
.82452 
.84420 
.86317 
.87528 
.8S096 
.89810 
.90890 
.92404 
.98348 
.04235 
.95065 
.95836 
.96545 
.97191 
.97491 
.97774 
.98041 
.66291 
.99570 
.99920 

-0.08630 
—.22400 
-.33200 
-.32383 
-.36922 
-.43392 
-.46284 
-.49015 
— . 51318 
-.52914 
-.53124 
-.63500 
-.63260 
-.52665 
-.61933 
-.50083 
-.48306 
-.45178 
-.44148 
—.41909 
40702 
-.38786 
-.86746 
-.35318 
-.33840 
-.32312 
-.30737 
-.28289 
-.26603 
-.24877 
-.23114 
-.21316 
—.19484 
-. 17623 
10082 
-.15735 
-. 14781 
-.13822 
-.06978 
-.02998 

20.00 

10.00 

5.00 
3.90 
3.10 
2.40 

2.00 
L60 
1.34 
1.18 
1.06 

.94 

.88 

.82 

.78 

.74 

.70 

.66 

.62 

.60 

.66 

.52 

.48 

.44 

.40 

.36 

.34 

.32 

.30 

.28 

.26 

.24 

.22 

.20 

.16 

.10 

.06 

.04 

0. 19603 
.39210 
.78431 

1.00533 

1.26502 

1.03399 

1. 90078 
2.45098 
2.92054 
3.32330 
3.69959 
4.17188 
4.45638 
4.78240 
5.02705 
5.29942 
5.00224 
5.94177 
6.32511 
6.53695 
7.00280 
7.54148 
8.16893 
8.91200 
9.80392 

10.8932 
11.6340 
12. 2549 
13.0719 
14.0050 
1&0630 

16.3399 
17.8253 
19.0078 
24.5098 
39.2167 
05.8595 
98.0392 

a 01042 
-.02790 
.13630 
.16896 
.19688 
.32318 
.44414 
.69013 
.68966 
.74934 
.79216 
.83260 
.86166 
.86975 
.88139 
.89258 
.90331 
.91866 
.92332 
.92801 
.93699 
.94543 
.95331 
.96062 
.96734 
.97345 
.97630 
.97897 
.98160 
.98386 
.98607 
.98812 
.99001 
.99173 
.99470 
.69793 
.99926 
.99967 

-a 06474 
-.18977 
-.33799 
-.36747 
-.44344 
-.63919 
-.66786 
-.65477 
-.61771 
-.48289 
-.46080 
— . 41378 
-.39360 
-.37209 
-.35722 
-.34188 
-.32610 
-.30988 
—.29328 
-.28480 
-.26760 
-.25004 
-.23214 
-.21393 
-.19543 
-.17686 
-.16718 
-. 15765 
-.14806 
13842 
-.12873 
-.11900 
-.10923 
-.09942 
-.07971 
-.04993 
-.02908 
-.02000 

20.00 

iaoo 

5.00 

3.10 
2.50 
LOO 
1.60 
L30 

1.10 
.94 
.84 
.76 
.68 
.62 
.68 
.64 
.60 
.48 
.44 
.42 
.40 
.38 
.36 
.34 
.32 
.30 
.28 
.25 
.24 
.22 
.20 
.18 
.16 
.06 
.06 
.04 

a 15625 
.31250 
.62500 
1.00806 

1.25000 

I. 64474 
L 95312 
2.40385 
2.84091 
3.32447 
8.72024 
4.11184 
4.59559 
5.04032 
5.38793 
5.78704 

6.25000 
0.51042 
7.10227 
7.44048 
7.81250 
8.22308 
8.68056 
9.1 OILS 
9.76562 

10.4167 

II. 1607 
12.0192 
13.0208 
14.2045 
15.6250 
17.3611 
19.6312 
39.0625 
52.0833 
78.1250 

I 

-0.07002 

-.13439 

-.32400 

-.49133 

-.57452 

—.00118 

-.57738 

-.62374 

-.47109 

-.41970 

-.38309 

-.36288 

-.32043 

-.29509 

-.27776 

-.26009 

-.24211 

-.23301 

-.21460 

-.20030 

-.19594 

-.18651 

-.17703 

-.16760 

-.16791 

-.14828 

-.13800 

-12SSS 

-.11912 

-.10932 

-.09949 

-.08963 

-.07974 

-.03996 

-.02938 

-.02000 

Af=2 


TU 10 

Af-S 

20.00 

10.00 

5.00 
2.70 
2.10 

1.00 
1.30 
L 10 

.60 

.80 

.74 

.70 

.64 

.68 

.54 

.50 

.46 

.42 

.40 

.38 

.36 

.34 

.30 

.28 

.26 

.24 

.22 

.20 

.18 

.16 

.14 

.10 

.06 

.04 

.02 

Si 

1 

1 

mm 

Kill 

IP nT!*r 

Bj 

pgij 

|g utW 

20.00 
10.00 
6.00 
4.80 
2.40 
LOO 
L 40 
L20 
L04 
.96 
.86 
.78 
.72 
.68 
.62 
.58 
.52 
.48 
.46 
.42 
.38 
.36 
.34 
.32 
.30 
.28 
.26 
.24 
.22 
.20 
.18 
.16 
.14 
.12 
.06 
.04 
.02 


m 

IB 

■ ' Atyir 
gin jffCT 

|J|B 

IS ’■*;**,? 


2a oo 
iaoo 

5.00 
4.40 
2.20 
1.80 
1.30 
1. 10 
.88 
.80 
.72 
.66 
.62 
.68 
.62 
.46 
.44 
.42 
.38 
.36 
.34 
.32 
.30 
.28 
.28 
.24 
.22 
.20 
.18 
.16 
.14 
.12 
.10 
.06 
.04 
.02 

0.10889 
.21978 
.43956 
.49950 
.99900 
1.22100 
1.69062 
1.99800 
. 2.49750 
2.74726 
3.05260 
3.33000 
3.54484 
3.78931 
4.22654 
4.77783 
4.99500 
5.232S6 
5.78369 
6.10501 
5.46412 
0.86813 
7.32601 
7.84929 
8.45309 
9. 15761 
9.99001 
10.9S90 
12.2100 
13.7363 
15.6986 
18.3150 
2L 9780 
35.6300 
54.9451 
109.890 

R 

IB 

B| yi'V/ 

-a 05304 
-.08576 
-.21422 
-.31881 
-.68957 
-.65969 
-.56336 
-.49023 
— .408S2 
-.37644 
-.34274 
-.81666 
-.29891 
-.28091 
-.26343 
-.22544 
-.21000 
-.20662 
— . 18742 
-. 17781 
16615 
-.16846 
-.14873 
-.13896 
-.12917 
11935 
-. 10950 
-.09062 
-.08972 
-.07981 
-.06987 
-.06992 
-.04995 
-.02998 
-.02000 
-.01000 

2a oo 
10.00 

5.00 
4.20 
2.10 
1.70 
L 20 

1.00 
.84 
.80 
.72 
.68 
.66 
.62 
.66 
.60 
.46 
.42 
.38 
.36 
.34 
.32 
.30 
.28 
.25 
.24 
.22 
.20 
.18 
.16 
.14 
.12 
.10 
.06 
.04 
.02 

a 10417 
.20833 
.41667 
.49603 
.99206 
1.22549 

I. 73611 
2.08333 
2. 48016 
2.60417 
2.89352 
3.06373 
3.15657 
3.38022 
8.72024 
4. 16667 
4.52899 
4.96032 
5.48246 
5.78704 
0.12745 
6.51042 
6.94444 
7.44048 
8.01282 
8. 68056 
9.46070 

la 4167 

II. 5741 
13.0208 
14.8810 
17.3611 
2a 8333 
34.7222 
52.0633 

104.167 

. 

BES 

■ 

■ 
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REPOET NO. 840 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TABLE II.— VALUES OF FUNCTIONS USED IN THE FLUTTER CALCULATIONS 
The expressions employed in the calculations of this table are: 


{- 2CWr+ t [* (a?) sfn =“S Jl (b) 008 “]} 
L '~ |“ 2 Uo)/+x [• Jt (s) ** “+B Jl (b) sin = ]} 


L\’ — ii+Ai 


U f -U~ Li+A* 


Al ~ jMuri [m B Ji (.v) 008 “~ Ji (jf) Sln 3 ] 
Ai ~ (sr) (|f) «»«] 


Afi'-Zi-Ai 


(ij— jBi 3 — j- CEt-f-ili) 


B '~ [ 3 Jl (s) cos3+ Af ; « (m) ws “+' 7 > (at) 8ln= ] 

a "^izfsMI Jl (I) (I) "»«+•» (I) “>=] 


L*' 

AT,' 

Aft' 

A£t' 


Aft'+i,' 

Afi'+At' 

Dn 

1>I 


82.4688 
136. 137 
218.902 


6.92663 
11.7828 
27. 4433 
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TABLE II— VALUES OF FUNCTIONS USED IN THE FLUTTER CALCULATIONS— Continued 




0.02700 
.10722 
.456 SO 
. 73408 
L 14404 
L9S049 
2.99507 
4.97860 
7.40271 
0.78096 
12.3445 
15.9586 
18.3504 
21.2901 
23.6401 
26.3831 
29.6115 
83.4462 
3S.04S8 
40. 7032 
46.8924 
54.5661 
64.2390 
76.6703 
03.0168 
116. 112 
129. 196 
146.004 
186.2S4 
191.064 
221.781 
260. 493 
310.237 
375.641 
5S7.618 
1506.20 
4186.04 
942a 14 


0.136SS 
.28916 
.51311 
.55892 
.59322 
.78890 
L 07377 
L 68663 
2.11529 
2.55920 
2.97716 
3.49574 
3.80467 
4.15586 
4.41806 
4.70685 
6.02866 
6. 38297 
5.78275 
6.00162 
6.48398 
7.03721 
7.07886 
S. 43288 
9.332S8 
10. 4274 
11.0689 
1L7891 
12.6035 
13.6324 
14.6021 
15.8479 
17. 8178 
19.0791 
23.9120 
38.3686 
64.0217 
96.0801 


a 02705 
.10313 
.45040 
.69894 
1.02716 
L 74282 
2.67381 
4.57135 
6-94229 
9.29616 
1L8325 
15.4270 
17.8097 
20.740S 
23.0853 
25.8230 
29.0463 
32.8763 
37.4742 
40.1266 
46.3117 


92.4224 
114. 614 
128.598 
145.404 
165.684 
19a 463 
221.179 
259.890 
809.634 
875.032 
687.005 
1505.59 
4185.51 
9418.46 



a 03143 
.08997 
.38710 
.63103 
L13S66 
2.22240 
3.44057 
5. 65117 
8.22275 
10.6933 
13. 3119 
16.9830 
19.4015 
22.3606 
24.7328 
27.4913 
3a 7346 
34.5835 
39. 1994 
4L 8604 
48.0620 
55.7473 
66. 4310 
77.8723 
94.2280 
116.332 
130.419 
147.231 
167. 516 
192.298 
223.018 
26L 733 
311.479 
376.885 
588.867 
1507.45 
4187.30 
9421.38 


-0.0064S 
-.02141 
-.07021 
-. 12311 
-.25428 
-.55617 
91682 
-1.59692 
-2.410S6 
-3.20292 
-4.04780 
-5.23712 
-6.02275 
-6.98717 
-7.75737 
-8. 6558S 
-9.71303 
-10.9077 
-12.4734 
-13.3416 
-15.3666 
-17.8745 
-21.0362 
-25.1030 
-30.4485 
-37.6635 
-42.2590 
-47.7516 
-54.3875 
-62.4490 
-72.53S8 
-85.2274 
-10L310 
-122.731 
-192.393 
-492.234 
-1395.70 
-2778.09 


1.57362 
L 92650 
2.31278 
2.62575 
2.93308 
3.31011 
3.65468 
3.92238 
4.22983 
4.58479 
4.78400 
5.23577 
5.49340 
5.77642 
0.08882 
6.43551 
6.82254 
7.25745 
7.74982 
8.31196 
8. 95997 
9.71631 
10.6072 
11.6767 
12.9830 
14.6148 
29.2800 
39.0498 
68.5353 


5.71955 
7.94412 
10.02S8 
12.3244 
15.4793 
18.6904 
2L407S 
24.7506 
28.9286 
31.4201 
37. 4613 
4L1494 
45.4047 
6a 3494 
56. 1411 
62.9848 
71.1517 
81.0058 
93.0463 
107.971 
126.781 
16a 951 
182.729 
226.679 
235.725 
1144.03 
2034.13 
4577.25 


.73691 
.81960 
.94318 
1.06737 
1.20910 
1.32773 
L 44693 
L 59609 
1.73469 
L 84380 
L 96975 
2.11660 
2.19951 
2.3SS60 
2.49701 
2.61648 
2.74870 
2.89612 
3.06104 
3.24691 
3. 45784 
3.69940 
3.97847 
4.30446 
4.09027 
5-15358 
5.72058 
6.43004 
12.S2SS 
17.0995 
25.6355 


1.09136 
1.46616 
1.82901 
2.22533 
2.54562 
2.85924 
3.24294 
8.59276 
3.86456 
4. 17518 
4.53388 
4.73500 
5.19061 
5.45019 
5.73518 
0.04956 
6. 39822 
6.78726 
7.22420 
7.71866 
8.28278 
8.93280 
9.69021 
10.5842 
1L 6558 
12.9641 
14.5980 
29.2706 
39.0425 
58.5345 


0.01474 

a 12493 

.05672 

.24911 

.23936 

.49241 

.57390 

.76842 

.89651 

.88199 

1.65534 

1.02910 

2.43919 

1. 13661 

3.86735 

1.29650 

5.65996 

1.45784 

7.77701 

L 64279 

9. 85747 

1.79822 

12.1500 

L 95482 

15.3021 

2.15130 

18.5113 

2.33422 

2L2270 

2.47842 

24.5692 

2.64505 

28.7461 

2.83956 

3L 2371 

2.94944 

37. 2778 

3.20019 

40.9650 

3.84409 

45.2197 

3.50264 

50.1641 

3.67846 

55.9556 

3.87415 

62.7989 

4.09350 

ra 965i 

4.34048 

80.8193 

4.62114 

92.8591 

4.94261 

107.784 

6.31420 

126.593 

5.74S28 

150.763 

6.26151 

182.540 

6.87952 

225.491 

7.63546 

285.537 

8.6S106 

1143.88 

37.1434 

2033.99 

22.6518 

4576.81 

34.8559 


a 01273 
.05987 
.18619 
.59724 
1.05363 
1.99355 
2.87470 
4.39070 
0. 14874 
8.40533 

10.5093 
12. 8191 
15.9871 
19.2071 
21. 9296 
25.2781 
29.4610 
3L9548 
38.0004 
41.6906 
45.9478 
6a 3944 
56.6379 
63.5383 
71. 7017 
8L 6573 

93.5093 
108. 525 
127.336 
161.603 
183.287 
226.238 
286.284 

1144.69 

2034.69 
4577.82 


—0 00294 
-.01098 
-.03415 
-. 11745 
-.21271 
-.42283 
-.62980 
-.99761 
-1.42724 
-L9S627 
-2.60902 
-3.03400 
-3.87400 
-4.67758 
— 5.35735 
-6.19330 
-7.23871 
-7.86162 
-0. 37212 
-10. 2946 
— 1L3577 
-12. 5950 
-14.0424 
-15.7650 
-17.7943 
-20.2582 
-23. 2722 
-27.0066 
-3L 7120 
-87. 7401 
-45. 7121 
-56. 4691 
-71.4066 
-288. 132 
-496. 917 
-1222.90 
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TABLE II.— VALUES OF FUNCTIONS USED IN THE FLUTTER CALCULATIONS— Coatinued 



























FLUTTER AND OSCILLATING AIR-FORCE CALCULATIONS IN TWO-DIMENSIONAL SUPERSONIC FLOW 
TABLE II.— VALUES OF FUNCTIONS USED IN THE FLUTTER CALCULATIONS — Concluded 





0.03287 
.00520 
.13262 
. 14740 
.29222 
.36313 
.51431 
.61290 
.77266 
.86222 
. 04928 
L 03736 
1. 10548 
L 18203 
1. 32129 
L 49554 
1.66413 
L 63023 
L 81305 
L 91440 
Z 02764 
2. 16500 
2.29929 
2.46410 
2.65435 
2.87619 
3.13830 
3.45278 
3.83708 
4. 31737 
4.93480 
5.76794 
6. 91021 
1L 3187 
17.2788 
34.5536 



0.00363 

0.04398 

0.00337 

.01453 i 

.08800 

.01412 

.05732 

.17594 

.04894 

.07360 

.19960 

.06314 

.29991 

.38813 

.31488 

.45309 

.47080 

.47689 

.88113 

.64818 

.91539 

L 28704 

.76184 

L 27497 

L 94260 

.94891 

198406 

2.35434 

L 04283 

Z 39693 

% 91092 

L 15771 

Z 05454 

3.46780 

1.26222 

a 61213 

3.93221 

L 34316 

3.97668 

4.49601 

143530 

4.54120 

5.59811 

L 60015 

5.64387 

7. 15913 

L80S09 

7.20542 

7.82656 

L 89003 

7.87302 

a 69161 

197978 

8.63822 

10. 4999 

Z 18767 

10.5469 

11.7012 

2. 80895 

117483 

13.1208 

2.44452 

13. 1670 

14. 8147 

2.59705 

14.8619 

16.8586 

2.76993 

16.9059 

19. 3558 

Z 96753 

19.4033 

22.4514 

a 19555 

22.4990 

26.3527 

3.46158 

26.4003 

31.3657 

3.77602 

814134 

37. 9667 | 

4.10336 

38.0044 

48.8648 

4.61461 

46.0126 

59. 3185 

5.19110 

59.3665 

77.4834 

5.93236 

77.6315 

105.471 

6. 92114 
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3.88152 
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3.54306 

615199 

4.52844 

61.6328 

4.25106 

SS- 5906 

5.43403 

88.6037 

7.08722 

246. 095 

9.05720 

246. 106 

ia6312 

553.716 

13.5864 

553.730 

21 2628 

2214. 88 

27. 0871 

2214.90 


1.98058 
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